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Abstract It is shown numerically that strange matter rings permit a con- 
tinuous transition to the extreme Kerr black hole. The multipoles as defined 
by Geroch and Hansen are studied and suggest a universal behaviour for 
bodies approaching the extreme Kerr solution parametrically. The appear- 
ance of a 'throat region', a distinctive feature of the extreme Kerr spacetime, 
is observed. With regard to stability, we verify for a large class of rings, that 
a particle sitting on the surface of the ring never has enough energy to escape 
to infinity along a geodesic. 
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1 Introduction 

In this paper, we consider in detail the parametric transition of a strange 
matter ring to a black hole. In [IIH], necessary and sufficient conditions for 
a quasi-stationary transition were presented and it was proved that an ex- 
treme Kerr black hole necessarily results. Using the analytic solution for the 
relativistic disc of dust [3] , a transition to a black hole was found explicitly 
[4] . Transitions have also been found numerically for rings with a variety of 
equations of state [S1[S] . Such parametric transitions to a black hole can be 
interpreted as a quasi-stationary collapse. 

The methods we use to study parametric transitions differ from those 
in the above cited papers, since we here concentrate on the behaviour of 
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multipole moments and on the appearance of a region of spacetime typical 
of metrics close to the extreme Kerr limit. These transitions are studied 
for strange matter, which is considered to be a form of matter that may 
be astrophysically relevant and has not been considered elsewhere for ring 
topologies. We restrict our attention here to ring configurations, since there 
is evidence suggesting that only rings and discs permit a transition to a 
black hole. We include in this paper a comparison with the corresponding 
transitions of rings governed by other equations of state. 

Section [5] is devoted to a brief description of the equation of state used 
here to model strange matter. In Sec. [3] we define multipole moments at 
infinity and follow their progression as they tend to those of the extreme 
Kerr black hole. The appearance of a "throat region" separating the "inner" 
from the "outer world" is discussed in Sec. ID In Sec.[5l we verify numerically 
that a particle resting on the ring's surface is always gravitationally bound, a 
condition, which can be considered to be a minimal requirement for stability. 
We close with a short summary in Sec. [H 

Throughout this paper, units are used in which the gravitational con- 
stant G and speed of light c are equal to one. 

2 Equation of state 

Strange matter is a fluid made of up (u), down (d) and strange (s) quarks. 
Our equation of state (eos) to characterize strange matter is the same one 
described by Gourgoulhon et al. 7J, who studied the properties of axially 
symmetric, stationary, spheroidal strange matter configurations. Based on 
the MIT bag model, we consider equal numbers of massless, non-interacting 
u,d,s quarks, confined to a given volume, i.e. enclosed in a "bag". The limits 
of the bag correspond in our case to the surface of the star, such that the 
star is entirely composed of strange matter. This model leads to a simple 
eos: 

e = ip + AB, (1) 

where e is the energy density, p the pressure and B the bag constant, char- 
acterizing the quark confinement. 

In the Newtonian limit, the pressure p is low and negligible in comparison 
to S, so the eos takes the form e — constant. Therefore, all the known 
Newtonian solutions for homogeneous bodies will be found in the Newtonian 
limit of the MIT bag model. Of course, a quark model of matter is not 
relevant in the Newtonian limit, but is taken as a limiting case of our eos. 
Also, like the homogeneous eos, but unlike polytropic models, the density 
of strange matter is discontinuous at the surface. 

3 Multipole Moments 

It has been shown in [1] that the extreme Kerr solution is the only black 
hole limit of rotating perfect fluid bodies in equilibrium. The extreme Kerr 
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black hole is characterized by the relation 

J = ztM^, (2) 

where M is the mass and J the angular momentum. To study quasi-stationary 
transitions that lead to black holes, we use bodies with a ring topology, 
since spheroidal bodies do not seem to have stationary sequences that lead 
to black holes [8 . For spheroidal bodies, a finite upper bound is observed 
for Zq, which is the relative redshift of zero angular momentum photons 
emitted from the surface of the body and observed at infinity. In contrast, 
the transition to a black hole occurs if and only if ^ oo [2] ■ We explore 
here such transitions with the concept of multipole moments. 

3. 1 The Metric and the Definition of Multipole Moments 

The line element for a stationary and axisymmetric spacetime containing a 
uniformly rotating fluid can be written in the form 

ds^ = e~^^[e^''{dp^ + dC^) + W^dip^] ~ e^^{adip + dtf, (3) 

where the functions e^'', e^*^, W and a depend only on p and C,. The equa- 
torial plane is given by C = and the axis of rotation by p = 0. 

To describe the surface of the ring, it is useful to introduce the metric 
potential V, 

where J7 is the angular velocity of the fluid with respect to infinity. The 
function V is constant along isobaric surfaces and the surface of the ring, 
defined to be the surface of vanishing pressure, can thus be denoted by 
V = Vq. The constant Vq is related to the relative redshift Zq via 

e-^« -l^Zo. (4) 

Consider now the vacuum region exterior to the mass distribution and 
extending to infinity. It is not possible to solve Einstein's equations in this 
region alone since the boundary conditions valid on the surface of the ring 
(and indeed the location of this surface) can only be found after solving the 
global problem. Considering the vacuum region will suffice for introducing 
the multipole moments however. In this region, there exists a conformal 
coordinate transformation z' — z'(z) {z' := p' 4- iQ' , z :— p + allowing 
one to choose p'{p, C) — W{p, C), which then leads to the metric 

ds^ = e-2^[e2fc'(dp'2 -f dC'^) + p'^dp^] - e2^(ad(^ + dt)\ (5) 

which we will use below to define the multipole moments. Note that the 
Cauchy-Riemann conditions for the transformation from ([3]) to (O imply 
W^pp +W,(^Q = 0, which is valid only in the vacuum domain. It follows from 
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axial symmetry that W — holds for p = {W = 0{p)). Along the axis of 
rotation, p = 0, one of the Cauchy-Riemann conditions then yields 

W 

lim — . (6) 

After solving numerically for the metric functions in Eq. an expansion 
of W/ p along the axis of rotation then allows us to find an expansion of 
in terms of C and vice versa. Thus, taking into account that p = <^ p' = 0, 
we are in a position to be able to write down the series expansion about the 
point C = +00 for the metric functions, which will be used in defining the 
multipole moments (see Eq. [T3| . 

Turning our attention back to Eq. (O, the Einstein equations governing 
a and e^^ can be rewritten using the single, complex Ernst equation 

(5i/)A/ = V/ • V/, (7) 

where / is the complex function / = e^^ + ?6, SR/ is the real part of / 
and A and V are respectively the Laplacian and the gradient operators in 
a three dimensional Euclidean space. Once a and U have been solved for, 
the metric function k' can be calculated via a line integral. Solutions of the 
Ernst equation lead to solutions of the Einstein equations and the metric 
potentials can be calculated from: 



P'e-'^b,c, (8) 



a<' = -pe b,p' (9) 
k'„, = p'[U,l, ^U,l -h,l )] (10) 

k',c' - 2p'[U,p, U,c' +^b,p, b,c' ]■ (11) 
From the Ernst potential /, one can define another potential ^: 

Taking the potential ^ on the positive part of the axis of rotation [p' — 0, 
C > 0), we can make a series expansion of it at infinity: 

00 

e(p'==0,C') = E7ST- (13) 

ri=0 ^ 

We assume reflectional symmetry about the equatorial plane in this paper 
for which it follows that m„ is real for even n and imaginary for odd n [HI 

The multipole moments P„ defined by Geroch [11] and Hansen [12] are 
algebraic combinations of the coefficients m„ and characterize the Ernst 
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potential uniquely. In this article, we consider the 7 first multipole moments 
of the infinite set given in [13] as 
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where Mjk = rrijnik — rrij^imk+i. Wc point out that Pq = Ad and Pi — iJ 
always hold. The multipoles P„ can then be normalized as follows: 

2/„ = t{-2tnr+^Pn. (15) 

For the Kerr black hole, the multipole moments are simply 

p^Kerr) ^ M{iJ/M)'', (16) 

where M and J are respectively the mass and the angular momentum of 
the black hole. Using the relation 

1 + 4122^2' ^^'^ 
where J7h is the angular velocity of the horizon, we then find 

2/^'-^H2/o)=Z/o(^)". (18) 

Through this normalization, all multipoles j/„ of the extreme Kerr black hole 
are equal to one, as can be seen by taking into account j/o = ^fi^iM = 1. 



3.2 Aiultipole Moments of Rings 

As Vb tends to — oo, we expect the multipole moments to become closer 
and closer to those of an extreme Kerr black hole. We tested this numer- 
ically by making use of a (slightly modified version of a) highly accurate 
computer program as described in [14j . This program was used for all the 
results presented in this paper. Since the n}^ multipole moment requires the 
calculation of n + 1 derivatives, which results in a loss of accuracy, we could 
not calculate arbitrarily many numerically. 

Figures[T]and[2]show the first seven multipole moments for homogeneous 
and strange matter rings where the ratio between the inner coordinate ra- 
dius Pi and the outer radius po (see Fig. [3]) is held constant at a value of 
Pi/ Po = 0.7. The left side of the plots corresponds to the Newtonian limit 
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Fig. 3 Example of a meridional cross-section of a strange matter ring. The ring 
in this example has the parameters Pi/po = 0.7 and eP^° = 0.1. 
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and the right side tends toward the black hole limit. As Vq — oo, the nor- 
malized multipoles all tend to one, demonstrating that this sequence indeed 
approaches the extreme Kerr solution. 

It is interesting to note, with respect to (or Zq), how slowly the 
exterior spacetime approaches that of a Kerr black hole. Consider, for ex- 
ample, the configuration from Fig. [2] with — 10~^. Whereas the value 
J/IVP = 1.00014 is very close to the limiting value of one reached in the 
extreme Kerr limit, the product 2QM — 0.9813 deviates rather significantly 
from it. This makes itself felt particularly for the higher multipole moments 
where powers of J7 are in play. The moment ?/4, for example, has reached 
only a value of yi^ ~ 0.91 for this configuration. 

To understand better the nature of the transition to the black hole, we 
compare the multipole moments of the above strange matter ring sequence 
with those of the Kerr solution. In Fig. 2] the y„ for n = 1 ... 6 are plotted 
against yo — 2f2Af for the strange matter ring sequence from above. A cor- 
responding picture for the sequence of Kerr solutions (see HH])) is displayed 
in Fig. m The clear similarity between these plots is emphasized in Fig. [6] 
where each y„ for the ring (solid line) and the Kerr solution (dotted line) is 
compared in a small figure over its whole range. The region very close to the 
extreme Kerr limit is then shown for j/i-j/s in detail. The graphs strongly 
suggest that the slopes 

^(yo = 1) (19) 

dyo 

are the same for the Kerr family and for the strange matter ring sequence 
discussed here. In fact, we found these slopes to be indepedent of the specific 
eos being used0 For the Kerr solutions, it follows from (fT8|) that 

^(yo = l) = n+l, (20) 
dyo 

which leads us to the conjecture that (|20p holds true for all sequences of 
rotating bodies that admit the transition to an extreme Kerr black hole. This 
conjecture provides a universal growth rate with which the y„ approach 
unity. It would be interesting to test this conjecture for the analytically 
known rigidly rotating disc of dust and work in this direction is ongoing. 

In Table [1] a comparison of the values of the first five moments ?/„ for a 
variety of configurations all with e^" = 10~^ is provided. The set of config- 
urations chosen includes rings with various different eos and various radius 
ratios and also includes the uniformly rotating disc of dust. A discussion of 
the multipoles of this last configuration as well as plots analogous to Fig. [2] 
can be found in [15j . Since all multipole moments tend to one in the limit 
Vq ^ — oo, these multipoles will provide almost no way of distinguishing 
between various configurations close to this limit. 



^ We checked this for ring sequences governed by homogeneous, polytropic and 
Chandrasekhar eos as well as for the rigidly rotating dust family. The Chan- 
drasekhar eos describes a completely degenerate, zero temperature, relativistic 
Fermi gas. 
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yi 


y2 


y-A 


?/4 


s.m. (n/ro 


= 0.6) 


0.982 


0.964 


0.947 


0.930 


0.913 


s.m. (ri/ro 


= 0.7) 


0.981 


0.963 


0.945 


0.928 


0.910 


s.m. (n/ro 


= 0.8) 


0.981 


0.962 


0.943 


0.925 


0.907 


hom. (n/ro 


= 0.7) 


0.981 


0.963 


0.945 


0.927 


0.910 


pol. (n/ro 


= 0.7) 


0.982 


0.965 


0.948 


0.931 


0.914 


rel. disc of dust 


0.984 


0.969 


0.953 


0.938 


0.924 



Table 1 The multipole moments j/n for various configurations, all with e ° = 
10^^. The abbreviation 's.m.' refers to a strange matter ring, 'hom.' to a homo- 
geneous ring, 'pol.' to a polytropic ring with the polytropic index n = 1 and 'rel. 
disc of dust' to the relativistic disc of dust. 
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In contrast, we present the multipole moments for configurations near 
the Newtonian limit (e~^° = LI) in Table [51 Here one can see that there 
is far more variation amongst the rings and that the disc of dust differs 
significantly from any of the rings. The values in the table also reflect the 
fact that strange matter has the same Newtonian limit as homogeneous 
matter. 
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2.22 


1.21 


1.04 


8.92 


7.44 


s.m. (n/ro = 0.7) 


2.09 


1.16 


1.02 


8.69 


7.43 


s.m. (n/ro = 0.8) 


1.92 


1.09 


0.978 


8.49 


7.46 


hom. (n/ro = 0.7) 


2.09 


1.16 


1.01 


8.68 


7.42 


pol. (n/ro = 0.7) 


2.14 


1.27 


1.11 


10.1 


8.75 


rel. disc of dust 


2.36 


1.73 


1.56 


12.5 


21.7 



Table 2 The multipole moments y„ for various configurations, all with e 
1.1 Zo = 0.1. The configurations are labelled as in Table[T] 



A further comparison of rings of various eos can be found in Fig. [7l 
Sequences of rings rotating at the mass-shedding Hmit, are plotted in a 
two-dimensional parameter space with 1 — e^" on the y-axis and pi/po on 
the cc-axis. The mass-shedding limit is reached when the path followed by a 
particle rotating at the outer edge of the ring becomes a geodesic. For a given 
eos, other ring configurations (i.e. not rotating at the mass-shedding limit) 
lie to the right of the corresponding curve. One can see that a transition to 
the extreme Kerr black hole is a generic feature of all rings considered here. 
The transition to spheroidal bodies exists for strange matter rings, but not 
for all eos. What is particularly striking is how close together the curves for 
strange and homogeneous rings remain right up to the black hole limit. This 
figure is a modified version of Fig. 1 of [6] . A discussion of the polytropic 
and Chandrasekhar eos can also be found in that paper. 



4 Throat geometry 

One of the most interesting features of bodies near the extreme Kerr black 
hole limit is the appearance of a throat geometry P^BlHj. In the limit, the 
throat separates the 'inner world', containing the ring, from the 'outer 
world'. The outer world is the asymptotically flat extreme Kerr spacetime, 
which is described by a single parameter and in which the horizon is located 
at the end of the infinitely long throat. On the other hand, the inner world is 
not asymptotically flat and is related to the outer world through its asymp- 
totic behaviour, which contains information about the one free parameter 
that uniquely describes the outer world. Any point in the outside world is 
inflnitely far away from any point in the inner world. For example, in the 
equatorial plane, one flnds that the radial proper distance 6 from the point 
p = to the point p is 

p 

V9pp\c=o '^P- (21) 

If one measures the proper distance from a point pi to a point p in the 
above equation, then it tends logarithmically to oo for the extreme Kerr 
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Fig. 7 The parameter space for rings with a variety of eos is considered in the 
(pi/po)-(l — e^") plane. Each eos is bounded on the left by the corresponding 
mass-shedding curve. 



black hole as pi ^ (the horizon in the coordinates used here is located at 

One way to represent the throat is to plot ^Jg^jM in the equatorial 
plane as a function of 5/M. Then, the throat appears as a plateau, i.e. 
a region appears in which the circumference of a circle of constant radius 
p — Pc tends toward a constant, independent of the radius pc- As the 
extreme Kerr black hole is approached, this region becomes infinitely long. 
Figure [5] shows the appearance of the throat for a sequence of strange matter 
rings with pi/po — 0.7 as the parameter tends to zero. Even in the first 
of these pictures (e^^" — 10~^), the highly relativistic nature of the ring is 
demonstrated by the fact that a small portion of the curve has a negative 
slope. That is, there exists a region of spacetime in which circles lying in the 
equatorial plane and centred about the origin have decreasing circumference 
with increasing radius. The last of these pictures is similar to Fig. 13 in [17] 
in which the 'inner world' is separated from the extreme Kerr solution by 
the infinite throat region. The proper distance between a point in what 
becomes the inner world (e.g. the outer edge of the ring p = po) and a point 
in what becomes the outer world (e.g. p = M) tends to infinity as e^^' — > 0. 
In a sense, we can say that the 'throat region' near the black hole limit 



^ It is thus to be expected that as Zq becomes large, the outer coordinate radius 
of the ring, as measured in units of mass for example, will become small. The shape 
of the ring itself does not change qualitatively as can be seen in Fig. 3 of [5] which 
depicts the shape of the ring along a sequence tending to the black hole limit. 
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Fig. 8 The function ^/g^ in the equatorial plane is plotted versus proper dis- 
tance, both normalized with respect to the mass M. In the throat region, ^/g^/M 
tends to the constant value 2. All four plots were made for a strange matter ring 
with a radius ratio p\/ po = 0.7 and with a value for e^^" as indicated. In the last 
plot, A gives the proper distance in the Kerr metric to the reference point p = M. 
Note that the proper distance between any point in the 'inner world' region and 
any point in the 'extreme Kerr' region tends to infinity as e^^° — > 0. 
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'swallows' the information as to what kind of configuration is sitting at the 
centre, as can be seen in Table [TJ 

The numerical 'inner world' solution was produced with a program that 
prescribes the asymptotic behaviour of the throat region (see |16j). Since the 
'asymptotically flat computer program' is capable of rendering rings with a 
relative redshift Zq well in excess of 100, the metric behaviour provided by 
this program can be used as initial input for the Newton- Raphson method of 
the 'inner world program' [T4|. The fact that such inital data converges to an 
inner world solution is in itself strong numerical evidence for the existence 
of a continuous transition to the black hole. 

5 Escape energy 

With the four- velocity and the Killing vector ^ ~ d/dt corresponding to 
stationarity, one can define the specific energy of a test particle with respect 
to infinity, i.e. the energy per unit mass, as 

E = -u'e., (22) 

which is a conserved quantity along any geodesic. 

With referring to the four-velocity of a particle resting on the ring's 
surface, E — 1 could be called the "escape energy" . If it is negative, then a 
sufficiently small perturbation will not suffice to induce the particle's escape 
to infinity on a geodesic, and it is referred to as gravitationally bound. In 
proving that Vq —^ —oo is a sufficient condition for reaching the Black Hole 
limit [2, use was made of the reasonable assumption that particles on the 
fluid's surface are gravitationally bound. One expects that this minimal 
requirement for stability will always be satisfied. We now proceed to verify 
this assumption for a large class of rings. 

Figures [9] and \W\ show the value of E along the surface of a variety of 
strange matter rings as it depends on radius. The radial parameter (p — 
Pi)/ {po — Pi) is chosen such that it runs from to 1 for every ring. In Fig. [H] 
curves are plotted for a constant value Pi/po — 0.7 and for varying Vq. We 
see that E tends to 1 in the Newtonian limit, which follows directly from 
Eq. (|22|) . Figure [TOl shows the behaviour of E for various values of Pi/po 
and constant Vq. Since conflgurations with small Pi/po do not exist when 
Vb becomes too negative (see Fig. [7]), we chose Vb to be in the Newtonian 
regime in order to be able to consider a wide range of values for the radius 
ratio. For every example considered in Figs [9] and llOi a maximal value at 
the outside edge of the ring in the equatorial plane is reached, just as one 
would expect. It is interesting to compare these results with the relativistic 
disc of dust for which E = 1 holds at the outer edge independent of the 
value of Zo [TB]. 

Focussing our attention now on the outer edge of the ring in the equato- 
rial plane, we see in Fig. [TT] how E depends on Vq for a sequence of strange 
matter rings with pi/po = 0.7. It is apparent that a maximum is reached 
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Fig. 9 The specific energy E versus (p - 
of strange matter rings with pi/po = 0.7 
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Fig. 10 The specific energy E versus (p — Pi)/(po — Pi) on the surface of a variety 
of strange matter rings near the Newtonian limit (e^^° = 0.9). 



in the Newtonian limit. For rings rotating at the mass-shedding limit, the 

value of E is also significantly smaller than one for small e^". The results 
for homogeneous rings are very similar and we can verify that E <1 holds 
(i.e. the escape energy is negative) for a large class of rings. 
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Fig. 11 The specific energy E versus e^" at the outer edge in the equatorial plane 
for strange matter rings with p-J po = 0.7. 



6 Summary 

It was shown numerically that a parametric transition exists from strange 
matter rings to the extreme Kerr black hole. Whereas it is known analyti- 
cally that the eos describing strange matter tends to that of a homogeneous 
body in the Newtonian limit, it was shown here that properties of configu- 
rations with these two eos remain similar right up to the limit Vq — *■ — cx3. 

Figs [Hand [2] suggest that the transition to the black hole is rather slow as 
g2Vo _^ This can be made more precise through the comparison with the 
Kerr solution, which leads us to conjecture that for every stationary rotating 
body permitting a transition to a black hole, the multipole moments y„ tend 
to one according to the formula 



dyn , . 

:^(yo = 1) 



1. 



It is expected that E < 1 always holds on the surface of a fluid body. 
Indeed, this inequality was used to prove that for rotating fluids the extreme 
Kerr black hole necessarily results if e^^" ^ [2J . We have verified that this 
inequality is correct for a large class of rings. 

As our knowledge of astrophysical collapse scenarios improves, it will be 
interesting to see how strong the connections can be to the quasi-stationary 
collapse considered here. 
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